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The persistent proximity of insulating and metallic phases, a puzzling characterestic of mangan- 
ites, is argued to arise from the self organization of the twofold degenerate Cg orbitals of Mn into 
localized Jahn-Teller(JT) polaronic levels and broad band states due to the large electron - JT 
phonon coupling present in them. We describe a new two band model with strong correlations 
and a dynamical mean-field theory calculation of equilibrium and transport properties. These ex- 
plain the insulator metal transition and colossal magnetoresistance quantitatively, as well as other 
consequences of two state coexistence. 
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Doped perovskite manganites Rei-xAxMnO^, where 
Re and A are rare earth and alkaline earth ions, show a 
rich variety of electronic, magnetic and structural phe- 
nomena and phases P, Q- Unusual effects, such as 
insulator-metal(IM) transitions both as a function of x 
and of temperature T over a wide region x ~ 0.5, or even 
as a consequence of isotope substitution (O^^ ^ 
colossal magneto-resistance (OMR) near Tjm and 'melt- 
ing' of the charge/orbitally ordered insulator {x ^ 0.5) 
into a metal in a relatively small(5-7 Tesla) magnetic 
field all suggest that metallic and insulating phases are 
always very close in free energy. This is also reflected 
in the ubiquitous coexistence (static or dynamic) of two 
'phases', one insulating with local lattice distortion and 
the other metallic without lattice distortion, with length 
scales varying from IQA° to WA° 0. 

These phenomena are due to the dynamics of the eg 
electrons of Mn constrained by three strong on-site inter- 
actions, namely electron lattice or Jahn- Teller (JT) cou- 
pling which splits the twofold eg orbital degeneracy, fer- 
romagnetic eg spin-t2g spin exchange or Hund's coupling 
Jh and eg electron repulsion U . The respective energies 
are Ejt — 0.5 eV, Jh — 2eV and U ~ 5 eV, compared 
to the eg electron intersite hopping t ~ 0.2 eV which 
sets the kinetic energy scale Q . Understanding their ob- 
served consequences is one of the major challenges in the 
physics of strongly interacting electrons. Earlier theo- 
retical attempts neglect one or more of these strong in- 
teractions and make further approximations; the predic- 
tions do not agree with many characteristics of mangan- 
ites. For example the ferromagnetic Curie transition in a 
number of manganites is from an insulator to a metal, for 
0.2 ~ a: ~ 0.5. However, a theory with just Hund's cou- 
pling, due to Furukawa |^, finds only a metallic phase, 
while Millis, MuUer and Shraiman '6!|, who additionally 
include the electron- JT phonon coupling g but treat the 
JT distortions as static displacements, obtain a metal- 
metal transition crossing over to an insulator-insulator 



transition as g increases. 

We propose and implement here a new approach which 
incorporates the crucial effects of all the three interac- 
tions and is based on a new idea, namely that of coex- 
isting JT polaronic and broad band Cg states, which we 
believe is the key to manganite physics. The idea and 
some of its consequences are described and calculations 
based on a new two band model are then outlined. 

We first discuss the effect of large JT coupling g on 
the initially twofold degenerate Cg orbitals at each lat- 
tice site. There is one superposition (labelled £) which, 
when singly occupied, leads to a polaronic state with lo- 
cal octahedral symmetry breaking Mn — O bond dis- 
tortion and energy —Ejt- Its intersite hopping is re- 
duced (for {Ejt /^i^o) ^ 1 where ujq is the JT phonon 
mode frequency) by the exponential Huang-Rhys 7J or 
phonon overlap factor 77 ~ exp {—EjT/2tiujQ). Since 
77 ~ exp{—5) for LaMnO^, and being a local quantity 
is not likely to change much on doping, the I polaron 
bandwidth 2D* ^ 2ztri ~ fcB(125i^) is small, and is ne- 
glected in much of this paper. At each site, there is nec- 
essarily another, orthogonal state (labelled h). Its energy 
is zero on the fraction x of sites where the t electrons 
are not present (hole-sites), and is C7 = {U -\- 2Ejt) on 
the i electron sites. The b electron hops adiabatically 
(ht^^ i-^Q^), in an annealed random medium of repul- 
sive £ sites and (for a homogeneous orbital liquid) forms a 
band whose re-normalizcd bandwidth 2D increases with 
X as well as with and H. For, the inhibition of b 
hopping due to large U is reduced when there are more 
hole-sites, and that due to large Jh is reduced when the 
t2g spin order is enhanced. We show in this paper that the 
simultaneous presence of £ polaron states and the b band, 
the change in the effective band width of the latter with 
X, T, H etc. and the consequent change in their relative 
occupation (leading to an insulator eg for Ejt > D and 
metal for Ejt < D ) are the basis of many phenomena 
seen in manganites. Our picture requires the persistence 
of local JT distortions well into the metallic regime, as 



2 



indeed seen by direct probes of instantaneous Mn — O 
bond length [8j. The (. polaron motion is anti-adiabatic 
{2D* <g; hujo) while the h electron motion is adiabatic 
{2D ^ ?iwo) ) the respective bandwidths are exponen- 
tially small and of order the bare value. By contrast, in 
earlier theoretical work where the JT distortion is 
treated as static, both £ and b electron dynamics become 
adiabatic so that the bands are of comparable width, 
leading to consequences which disagree with experiment. 

The system of correlated ^ and b states described above 
can be modelled by the Hamiltonian 

Ha = {-EjT - E 4^.. - A* E (1) 

i,(7 i(7 
{ij),a id 

In Eq.(l), the I polaron has energy {—Ejt), the b elec- 
trons hop between nearest neighbour sites with an effec- 
tive amplitude i , and U is the effective repulsion between 
£ polarons and b electrons of the same spin at a particu- 
lar site |lO|. The common chemical potential is /i. Hm is 
the magnetic part involving the eg spins si and t2g core 
spins Si, namely 

H„i = —Jh E Si ■ Si ~ Jf E ' ~ E ' -^(2) 

i {ij) i 

In addition to the familiar Hund's coupling Jh and Zee- 
man coupling to an external field H (neglecting the rel- 
atively smaller term fissl ■ H) we have included in Hm 
a new ferromagnetic nearest neighbour exchange Jp be- 
tween the t2g core spins. 

The new term arises from a virtual double-exchange 
process in which an £ electron at a site i hops 
quickly(adiabatically) to a nearest neighbour site j and 
back. In the limit of large U and Jh , and approxi- 
mating the t2g spins as classical, i.e., Si = SCli where 
rii are unit vectors, the energy shift to second order in 
tis {P/2Ejt) ^{Cli ■ Clj + 1) {nu{l-nj) + ni.j{l-n,)). 
The energy denominator 2Ejt arises from the unre- 
laxed JT distorted intermediate state, the ■^{^i ■ % + 1) 
factor from large Jh, and the occupancy dependent 
terms from large U. We have further approximated 
this by the simple exchange term in Eq.(2), whence 
Jf {P/{2EjtS^)) x{1 - x). 

The model Hamiltonian Hib (Eq.(l)) can be motivated 
starting from a lattice model with two Cg orbitals per 
sited 0] Let al^ create an electron in orbital a at 
site i, { with |a = 1 >= \x^ — > and \a = 2 >= 
3^2 _ ^2 say), locally coupled to JT lattice modes 
{Qxi, Qzi) = {Qi, di) by the term i?}^, = 5aL^a/3a»/3 • Qi 
where f are the Pauli matrices. The eigenvalues of H'jrp 
are igQi , the corresponding electron creation operators 
being labelled 6^ and £'^ . The effective lattice potential 
energies for single electron occupation of these states are 



{K/2)Ql ± gQi where K is the phonon force constant. 
The lower (£) branch has a minimum at a static lattice 
distortion Qo = 9/K and energy —Ejt ~ —{g'^/2K). 
£\ in Eq.(l) are polaron creation operators given by rji£\ 
where rji { ^ exp{iQoPi/h) where Pi is the radial mo- 
mentum conjugate to Qi ) are unitary Lang-Firsov 
like displacement transformations. The exponential re- 
duction T] in the hopping amplitude of £ arises essen- 
tially from their ground state average. (We neglect the 
small fluctuations with respect to this average). No lat- 
tice distortion or reduction in hopping is associated with 
the higher {b) branch. The small leading intersite term 

H'a - {4bj + b\^j + /ic) describing ^ - 6 hy- 

bridization is not included in Eq.(l); some consequences 
are discussed later. Thus for {Ejt/Hcuo) » 1, and ne- 
glecting H'lf^, we are led to Eq (1). We have made further 
approximations which are realistic for a homogenous or- 
bital liquid, by statistically averaging quantities which de- 
pend on the orbital admixture angle {Oi}, eg the intersite 
hopping and the nearest neighbour Mn — O bond correla- 
tions (short-long), so that t and bj are to be regarded as 

averaged over {0i}. This is reasonable for 0.2 ^ a; ~ 0.5 
in most manganites, but poor for those other values of 
X for which one has strong orbital correlations or long 
range orbital order. Finally, in this regime, the chemical 
potential fi is chosen such that 

(nu) + (ubi) = fii + fib = {1 - x) (3) 

We have calculated the equilibrium and transport 
properties of (Eq.(l)) using the dynamical mean field 
theory (DMFT) jl^] which is exact in d = cjo, and is ac- 
curate in three dimensions if spatial correlations are not 
crucial. Ha, being a generalized form of the Falicov- 
Kimball mo del [l^ with additional spin interactions H„i 
(Eq.(2)) and constraint Eq.(3), is exactly soluble in this 
limit [13 • We treat the t2g core spins as classical as 
stated earlier, and their interaction in the Curie- Weiss 
mean field approximation, and work in the large Jh limit 
(i.e. SJn/t > 1). The local self energy Sii(w+) of the b 
electrons and the mean magnetization m = (fii) are de- 
termined self-consistently. From these we calculate the 
spectral density or density of states (DOS) P6(w) of the 
b electrons which determines their occupation relative to 
that of the £ state, the current-current correlation func- 
tion relevant for the Kubo formula for the electrical con- 
ductivity, etc[3- Most of our calculations use a semicir- 
cular bare DOS with a bare bandwidth 2Do0|. 

Our results for Pb{ijj) for different values of x and T 
are shown in Fig.l. Since the effective bandwidth 2D of 
the b electrons decreases significantly as x decreases (for 
any sizeable U), the b band bottom is above the £ level 
for small x e.g x = 0.1 (Fig. la), and the low tempera- 
ture state is an insulator, ferromagnetic because of Jp- 
(We ignore the anti-ferromagnetic super-exchange as it 
is dominated by the much larger Jp for a; ~ .1 ). On 
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FIG. 1: Spectral density or 6-DOS Pb{^) for various values 
of doping X and temperature T. The effective £ polaron level 
is marked as a vertical line. Parameters chosen are Ejt = 
-0.5 eV, Do = 1.2 eV , tJ ^ 5.0 eV Jf = 2.22, meV. (a) 
X = 0.1, T = 0, /I = —Ejt; ferromagnetic insulator, (b) 
X = 0.3, T = 180/C(< Tc = 240K); ferromagnetic metal, 
(c) X — 0.3, T — 350K ; paramagnetic insulator. Full(dotted) 
lines correspond to up(down) spin DOS. Occupied band states 
are shown shaded. 
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FIG. 2: The resistivity p(T) oi Lai-^^Ca^MnOz {x = 0.33) as 
a function of temperature T. Calculated (full lines for realistic 
anisotropic hopping, dotted lines for a semicircular DOS) and 
experimental (circles, from ref.0|) results for H = and 
ITesla are shown. Parameters, chosen to fit Tc and p{Tc), 
are Ejt = -0.5 (SDOS),-0.6 eV(RDOS), 2D„ = 2.4 eV 
(SDOS), 2.44 eV (RDOS), U = 5eV, and Jf = 2.2,7 meV . 



increasing x, D increases as well, and beyond a critical 
Xc for which D — Ejt, the low temperature state is a 
ferromagnetic metal, as in Fig. lb. We note that while 
both b and I states are occupied in the metal, most elec- 
trons are in the latter, polaronic state since ni = 0.66 
and fib = 0.04. Tc is largely due to Jf', for, fib being very 
small, so is the b electron contribution to Tc via conven- 
tional double exchange. On increasing T at this x, the 
spins disorder, reducing the effective b electron hop- 
ping or D. Fig.lc shows the DOS at T = 350ii:(> TJ- 
It is that of a small gap semiconductor. The metal insu- 
lator transition occurs very near Tc because of the strong 
feedback between D and the magnetization m. The car- 
riers in the paramagnetic state are b electrons excited 
across the relatively small effective band gap, not ther- 
mally unbound small polarons. The mobile carrier frac- 
tion fib is very small in the metallic phase (fib — 0.06 
for parameters corresponding to LaQ/jSro ^MnOs) and 
decreases rapidly 0| to a minimum at Tc, whereas the 
electron density is 0.7 per site . This is exactly the hith- 
erto unexplained inference for rip f f from the observed 
Drude weight in this compound [iGj. 

Our theory also describes transport properties fairly 
well. Fig. 2 shows the electrical resistivity p{T) versus 
T for model parameters chosen to fit Tc and p(Tc) for 
LaosjCao.ssMnOs 0. The results for semicircular 
DOS and tight binding DOS are nearly the same. We see 
clearly the sharp paramagnetic insulator to ferromagnetic 
metal transition, the former having a calculated effective 
electrical gap ^sJJ. of 34 meV while the experimental 
value is 48 meV [1^ . The resistivity falls dramatically to 
about 2 m^lcm just below Tc and does not decrease much 
further thereafter, in contrast to the observed residual re- 



sistivity values which can be as small as 50 fiflcm . This 
is a consequence of our neglect of inter-site i coherence 
and is discussed later. The resistivity in a field of 7 Tesla 
is also shown in Fig. 2. The colossal decrease(CMR) is 
apparent [l^. 

The properties of manganites vary strongly and char- 
acteristically with the ionic species Re and A. This 
is natural in our model, since the relative balance be- 
tween metal and insulator is strongly affected by small 
changes in Ejt and Do and the carrier density depends 
exponentially on [D — Ejt). For example, keeping Ejt 
fixed (for a given x = 0.3 say), and varying 2Do, we 
find the variations in physical properties shown in Fig. 3. 
Tc, being proportional to Jp, increases roughly as . 
The Curie transition changes from insulator-insulator to 
insulator-metal to metal-metal, as is seen experimentally 
in the sequence Pro.jCaa.zMnOs , Lao. 7 C ao. 3 MnO 3 and 
Lao.rSro.sMnOs . Do is believed to increase 

in this order. We also exhibit the enormous vari- 
ation in the calculated fractional magneto-resistance 
{AR{H)/R{H)}j,^j.^ for H = 7 Tesla as a function of 
Tc; the agreement with experiment |20| is very good. 

Thus, in the regime 0.2 < x < 0.4, a wide range of 
physical properties for a number of doped manganites can 
be understood physically and described quantitatively by 
our theory jl3|- The ratio {Ejt /Do) broadly determines 
the systematics and Do the basic energy scale {Jh and 
U being close to 00 for our purposes). 

The intersite hopping of the £ polarons arising from 
term fl^j, with a characteristic temperature scale of 
T* ~ {zTji/ks) — 125K, has many important conse- 
quences. Below T* the £ states can form a coherent band 
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FIG. 3: Physical properties as a function of the d electron 
bandwidth 2Do, for fixed Ejt = 0.5 eV. Full line: ferromag- 
netic Tc vs Do- The resistive transition is from an insulator 
to an insulator(PI-FI) for small Do, from an insulator to a 
metal(PI-FM) for intermediate Do, and from a metal to a 
metal(PM-FM) for large Do, as indicated . Broken line: cal- 
culated fractional magnetoresistance {AR{H) / R(H))j, vs Tc 
for H — 7 Tesla . Experimental points are from |2Q]. 



due to hybridization with b states, whence the JT dis- 
tortion can become dynamic (time scale ~ h/ksT*), and 
also smaller self-consistently. The 6 electron scattering 
and the consequent resistivity then vanish at T = 0, and 
are nonzero only if static disorder is present. This can 
lead to a metallic state with a small residual resistivity 
or to an Anderson localized insulating state depending 
on the amount of disorder. The strong dependence of Tc 
on isotope mass A/q is also a consequence of £ polaron 
hopping, since its direct double exchange contribution to 
Tc will be proportional to i{l — x)exp{—EjT/'2fiujo) and 
ujo = y/K/M f). The isotope effect thus estimated is of 
the right size |21|. 

Spatial correlations and inhomogeneities can be inves- 
tigated in an extended version of the model in Eq.(l) 
by including the orbital angles 9i as additional degrees 
of freedom. The intersite hopping tij now depends on Oi 
and 6j. There are a number of anharmonic, steric, elastic 
terms, eg., nearest neighbour long-short Mn — O bond 
correlations, coupling between JT modes and strain, etc., 
which depend on these angles as well. These, and other 
factors such as the variation of Ejt with local ion size, 
can be included in Eq.(l), and questions such as short 
range correlation, orbital order, long range order in 
{flu ± flu) (total charge/relative charge order), disorder 
induced 'phase' separation, explored by going beyond the 
self-consistent single site DMFT used here. It would also 
be interesting to explore the idea of exponentially sepa- 
rated time scales and adiabatic - nonadiabatic crossover 



in other oxides and organic molecules with orbital degen- 
eracy and strong local symmetry breaking Jahn Teller 
coupling. 
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